arXiv: 1503.05149v3 [physics.flu-dyn] 4 May 2015 


On Exact Solutions of the Navier-Stokes Equations for 
Uni-directional Flows 

By F. Lam 

In the present note, we show that the uni-directional flows in a rectangular chan¬ 
nel and in a circular pipe are exact spatio-temporal solutions of the Navier-Stokes 
equations over a short time interval. We assert that the classical plane Poiseuille- 
Couette flow and Hagen-Poiseuille flow are time-independent approximations of the 
exact solutions if an appropriate initial velocity distribution at starting location is 
specified. Conceptually, there do not exist absolute steady flows starting from un¬ 
specified initial data. The classic experimental measurements by Poiseuille can be 
explained in terms of the evolutional solutions. In particular, the pipe flow does 
not have a time-independent characteristic velocity. The orthodox notion that the 
parabolic profile exists for arbitrary Reynolds numbers is unwarranted. 
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1. Introduction 

In the Eulerian description of the motion of an incompressible, homogeneous New¬ 
tonian fluid, the momentum and the continuity equations for fluid dynamics are 

-I- (u.V)u = V.u = 0, (1.1) 

where the velocity vector u = u(x, t) has the components {u,v,w), the scalar 
quantity p = p(x, t) is the pressure, the space variable is denoted by x = (x, y, z), 
and A is the Laplacian. The density and the viscosity of the fluid are denoted by p 
and pL respectively. The kinematic viscosity is v = p/p. These equations are known 
as the Navier-Stokes equations. They are derived on the basis of the continuum 
hypothesis. 

We seek the solution of dm as an initial-boundary value problem in the space- 
time domain denoted by Hx [0, t > 0], where H has a smooth impermeable boundary 
dD,. The initial condition is given by 

u(x, t=0) = Uo(x) X G £1, (1.2) 

and the no-slip boundary condition is 

u(x, t) = 0 X G cl£l. (1.3) 

The equations of motion (ED are often re-written in dimensionless form 

-f (u*.V)u* =-^Au* - Vp*, V.u*=0, (1.4) 

Re 

where Re denotes the Reynolds number with reference to a characteristic length 
scale and a characteristic velocity. A prevailing concept is that there exist some 
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basic flows, such as plane Poiseuille flow and the Hagen-Poiseuille flow in a circular 
pipe, which are “exact” solutions of the Navier-Stokes equations, though these basic 
flows are in fact the solutions of simplified equations of motion. Nevertheless, one 
important implication is that the basic flows exist at all Reynolds numbers. In 
fact, the stability of a flow depends essentially on this idea which quantifies the 
growth or decay of disturbances as the Reynolds number varies, as a parameter. 
In the present note, we show that approximating or simplifying the equations of 
motion without due justifications can be very misleading, as the solutions of the 
approximated equations are inconsistent in the spatio-temporal evolution of fluid 
motion. 


2. Some particular solutions 

Let us consider plane Poiseuille flow in a rectangular channel with half width b. As 
in many textbooks on fluid mechanics, we make the following assumptions for the 


flow: 

V = w = 0, du/dt = 0, u = u{y), and p = p{x). (2-1) 

The momentum equation becomes an ordinary differential equation 

d^u/dy^ = {dp/dx)/pL, u(±&) = 0. (2.2) 

The time-independent steady solution of the boundary value problem is given by 

Alternatively, we seek Navier-Stokes solutions by assuming 

v = w = 0, u = u{t), and dp/dx = dp/dx{t). (2.4) 

Now the flow evolution is an initial-value problem: 

du/dt = —{dp/dx)/p, u{t = 0) = uq (2.5) 

Its solution defines a class of space-independent unsteady flows, 

= (|^)(r)dr, (2.6) 

which complements solution (12.31) by space-time symmetry. 


Both solutions (12.311 and (12.611 may be regarded as exact particular solutions of 
the Navier-Stokes equations, only in a restrictive sense. Moreover, plane Poiseuille 
flow is considered to exist at arbitrary Reynolds number Ret = Ub/v, where U is 
the velocity at the channel centre and > 0. Similarly, we claim that the unsteady 
solution (1^ is valid for all unit Reynolds number Re =\/v. 

Strictly speaking, the correctness of the profile (IQ) must be tested experi¬ 
mentally by intelligent beings squatting absolutely on space while supernaturals 
crawling in time may be able to verify the flow (EH)! 

For a straight circular pipe of radius A, we have the analogous expression to 

EH: 
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This is known as the Hagen-Poiseuille flow. 

In fact, we have no a priori justifications to detach space from time or vice versa, 
as implied in the simplification or (1^ . Every fluid motion develops in space 
as well as in time. Since we have separately considered the initial-boundary data 
and (El, neither (lO) nor (12.61) can be regarded as correct descriptions of 
the evolution. Practically, the velocity profile (IQ) is accepted as “correct” because 
there are experimental verifications (Hagen 1839; Poiseuille 1840) even though these 
experiments ought to have been carried out locally. Within the experimental errors, 
the parabolic profile agreed with the measurements of ensemble averages. The ex¬ 
periment repeats and the data sampling must have been undertaken over time, 
though the time-dependence may not be shown in an explicit formulation. 

We recall that the first theoretical explanation of (EB in circular tubes was 
not done by solving equations but was derived by Poiseuille by considering the 
momentum balance (see, for example, §331 of Lamb 1975) over a fully developed 
flow or a transient with no reference to its initial condition. The question is, to 
what extent, can we ignore the temporal effects? 

3. Flow evolution over space-time 

Instead of EB, we postulate 

Vt € [ts, ts + e]) e > 0, v = w = 0, u = u{y,t), and p = p{x,t), (3-1) 

where tg stands for starting time, and e denotes a short interval. The start refers 
to either the initialisation of a motion from rest or a subsequent moment when the 
assumed uni-directional flow has been established at least over a part of the appa¬ 
ratus. The second option is more relevant to experimental measurements because 
the flow at the inlet may be specified as initial data. With these understandings, 
we set tg = 0. 


(o) Plane Poiseuille flow 

We consider the following initial-boundary value problem: 

du/dt — vd'^u/dy^ = q{x,t), u{y,0)=uo{y), m( 0) = u(l) = 0. (3.2) 

where q{x,t)= —{dpldx)lp is the scaled pressure gradient (6 = 1). The solution of 
this diffusion problem is well-known, (see, for example, Carslaw & Jaeger 1947), 

u{y,t)= uo{a)H{y,a,t)da + / q{x,T)H{y,a,t-T)dadT, (3.3) 

Jo Jo Jo 

where Green’s function is given by 


H{y,<T,t) = 2 E exp ( — sin (mry^ sin (mra ), 


(3.4) 


n—1 


for all time t € [0, e], and in general e = e{uo,q). For constant pressure gradient 
q{x) = qo, the second integral in (13.31) is found to be 


4go 


sin(n7rj/) 


odd n 


^1 — exp ( — n^TT^vtj'^ . 


(3.5) 
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Consider the particular initial profile at a location closed to the inlet 

uo{v) = ay{l-y), (3.6) 

where a is a given constant, its contribution to the solution is given by 

H exp ( - n ^ vt). (3.7) 

odd n 

We make the following assertions: 

• In view of (13.51) and (1X71) . the subsequent flow development of initial data 
(13.6|) is steady over a short time interval if 

a = qo/{2v). 


The solution can be expressed as 


u{y,t) 


u{y) = 




odd n 


sin(n7r?/) 

n^TT^ 


1 

2fj, 





(cf. (12.31) 1. by virtue of the Fourier series expansion for functions y and y'^ 
over 0 < 2 / < 1. 


• For a ^ —{dp/dx)/{2y.) in (13.61) . the flow is time-dependent. However, the 
following approximation is valid. Within a short time f —>■ 0"*" or for small 
viscosity v such that the exponential decay has virtually no influence. Mea¬ 
surements made at a small distance downstream must be very close to the 
initial data, as implied in (13.71) . The small viscosity requirement is relevant to 
many experiments using air (v 0(10“^) in SI units, m^s“^) and pure water 
{v ^ 0(10“®)) under the conditions of standard temperature and pressure. 


• For an arbitrary initial velocity uq of a polynomial satisfying the boundary 
conditions, the flow is unsteady in general. If we ignore the contribution from 
the pressure gradient over the short-time interval e, the flow downstream of 
the inlet resembles the initial data, in the limit of i/t 0, according to Fourier 
series expansions for simple functions. 


(&) Plane Couette flow 

If we consider the following modihcation of the boundary conditions in dSJ): 

u{y = 0) = H(t), u{y = 1) = H(t), (3.8) 

we obtain the solution for plane Couette flow: 


i'{y,t)= / uo{cT)H{y,a,t)da + / / q{x,T)H{y,a,t-T)dadT 

Jo Jo Jo 

+ [ {A- B){t) h{y,t-T)dT, 

Jo 


(3.9) 
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where 

CXD 

h{y,t) = 2 E( — I)”n7r exp ( — n’^n^h't) sin (mry). (3.10) 

n—1 

If the boundary conditions are time-independent, the last integral can be evaluated 
as 

(3.11) 

n=l 

Let us consider several particular cases of our general solution: 


• An impulsively started motion mo= 0, with zero pressure gradient and bound¬ 
ary conditions A = 0, B/iy = -|-1. 

In view of (EH) and (I3.11|) . the time-independent solution of this motion is 
plane Couette flow, 

uiyG)'-uiy) = 2'£i-ir+^"-^^^^^=y, 0<y<l. (3.12) 

rnr 


Evidently, the viscosity plays a key role in determining how soon the temporal 
effects may be neglected. 

• A motion driven by the upper boundary B = Q sm(u}t)/2 (A=0) with initial 
data Mo=0 and zero pressure gradient. 

The solution is given by 


n—1 


(mr) sm(mTy) 
up' + 




(3.13) 


where 


x(a;) = {PtP'v) sin(a;t) — uj cos{ujt) + uj exp ( — PPvt). 

In practice, this result may only be realised within a short period of time 
after the sinusoidal force has been applied since the non-linearity cannot be 
ignored, especially for large Q. 

• A motion driven by A = asin(a;o,t)/2 and B = /3sin(w^t)/2 with zero initial 
data and zero pressure gradient. 

In view of the last solution, we find 

OO 

u{y,t) = y](-l)”+^ {mr) Sin(n7ry) 

n—1 

The solution can be considered as a model for simulating the effect of wall 
vibrations in experiments if we choose small amplitudes a and jS. 

More complicated flows can be generated and analysed if we include non-zero 
initial data and time-dependent pressure gradients. 


«X(^a) _ /0 143 

oj'i + PPP uji + PPP )' 
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4. Hagen-Poiseuille flow 

In a cylindrical co-ordinates system {r,9,z), consider (u,p) = {u,p){r,9, z,t). The 
velocity components in the radial and circumferential directions are hypothesised 
as identically zero over time t S [0, e]. The axial component is treated as a function 
of the independent variables r and t, 

w = w{r,t), p = p{z,t). 


By analogy, the differential equation for Hagen-Poiseuille flow in a circular pipe of 
radius unity becomes (see, for example, Batchelor 1973) 


dw ldw\ 

dt V dr"^ r dr ) 
w{r,0) = wo{r), 


w{l)=0, 


(4.1) 


where we seek bounded solution, |r/;(r)|<oo, Vr€[0,1), and ^ = — {dp/dz)/p denotes 
the axial pressure gradient. The solution of the initial-boundary value problem can 
be expressed in terms of the Bessel functions of orders zero and one. 


w{r, t) 


wo{s)G{r, s, t)ds + v 


■i/'(^, T)G'(r, s, t—r)dsdT, (4-2) 


Jo 

where Green’s function is given by 


0 JO 


r^l 4.\ r> / \2 >^0 do (A„s) 

G(r, s,t) = 2^ exp ( - s ' \ - 

n=l "'l vAnj 


where the constants, A„,n = 1,2, 
first few of them are 


are the positive zeros of Jo(A„) 


(4.3) 
0. The 


Ai = 2.40482556, 
As = 14.93091771, 
Ag = 27.49347913, 


Aa = 5.52007811, 
Ae = 18.07106397, 
Aio = 30.63460647, 


As = 8.65372791, 
At = 21.21163663, 
All = 33.77582021, 


A 4 = 11.79153444, 
As = 24.35247153, 
Ai 2 = 36.91709835. 


The other zeros may be calculated by making use of the asymptotic formula 


A„ =/3-I-1/(8/3), /3 = (n - l/ 4 ) 7 r. 


In particular, none of these zeros coincides with any zero of Ji{an)- Bessel functions 
J„(x) are entire functions of the argument x (Watson 1944), 


Explicitly, 


We notice that 




Mx) = i-- + 




22 22 42 22 42 62 


-h 


max|Jo| = l at a; = 0, max |Ji | w 0.5819 at a; « 1.8412. 
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As X tends to 0, Ji{x) tends to 0, and Ji(Ai) ss 0.5191. 

Evidently, the flow evolution reaches a steady state in the interval e if i/t —0+ 
when exp(—~ 0(1) for some n < 10 (say). The smallness assumption for i/t 
may be justified for air and water. Practically, the uni-directional flow can only be 
observed during an interval t ~ 0(1) measured in seconds under normal laboratory 
conditions. Thus the temporal effects due to the exponential decay reduces to the 
extent when any experimental measurements of the flow become time-insensitive. 

In view of the identity, 

/■i — 1V 

/ s(l - S^) —Vo(A„s)ds = - 

the integral involving the Bessel functions in Green’s function may be evaluated, 
particularly for the initial condition having a “parabolic” distribution. 

Assuming a time-independent pressure gradient, the second term in 63) con¬ 
tributes nothing as vt ^ Q. Thus the solution due to the initial data, 

wo{r) = I3{1 — (/? = const.) 


is given by 


w{r, t) 


w{r) = 4/? ^ 


n—1 


d2(A„) 

XI Jl{\u) 


JoiXn'^ ) 


a ST (_ _^2 I An 

"+16 

/3(1 -r^) -f 0(r'‘), 


r^- 


(4.4) 


where we have performed the necessary calculations. The leading constant term 
converges quickly to 1 for m > 7. The coefficient for is defined by an alternating 
series. We first choose m = 3,7, 15, 31 and extrapolate these results to m —)• oo to 
get 1.0620. We then repeat the process using m = 4,8,16, 32 to get 0.9405. By the 
minimal computational effort, we verify that the coefficient does equal to unity by 
considering the mean value (1.0013). 

In fact, if the initial velocity has a more general profile. 


k 

wo{r) =, (m > 0, finite fc), (4.5) 

i—m 


the time-insensitive solution can be computed in a similar manner. Consequently, 
we may contemplate that Poiseuille’s measurements ought to be the transient of 
well-controlled initial parabolic profiles, which can be generated and maintained by 
a fixed pressure gradient during the tests. 

Szymahski (1932) considered flow development from an impulsively started mo¬ 
tion from rest. He derived 


w{r, t) 


4/r dz 



CO 




Jo{Xk^ ) 

Afe'Ai(Afc) 


exp ( - Xlvt) 


(4.6) 
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where it is assumed that the pressure gradient is zero for t < 0, and approaches to 
a constant as t > 0. The Hagen-Poseuille flow is recovered as t oo. In reality, 
solution (ITBI) is misleading as it implies that the approximation dUD holds for 
arbitrarily long time. For air and water, the duration for the unsteady effect to 
disappear altogether must be in the order of t ~ 0(10®) seconds! On the contrary, 
high-Reynolds-number streamlined flows in pipes have only been observed in brief 
periods of time and over finite space. (Reynolds 1883; Pfenninger 1961; Mullin 2011), 
especially for fluids of small viscosity. 

One possible way out of the dilemma (14.61) is to adopt our idea of the limit 
i^t —!> 0+. The first term on the right is in fact the second term in (14.21) which is 
showed to be negligible in the short time limit. In practice, it is very difficult to 
apply a constant pressure gradient right from the beginning of the motion. The 
initial flow will be highly unsteady and its life span is likely too short to determine 
suitable modifications for subsequent repeats (so as to achieve a constant pressure 
gradient). However, it is relatively easy to produce established or full-developed 
flows, at least over a portion of the pipe, driven by an adjustable pressure gradient. 
We may continuously operate our apparatus to maintain a constant flow rate for 
sufficiently long time. Numerous experiences show that it is crucial to control the 
flow conditions near the pipe inlet in order to achieve a parabolic profile over useful 
distance downstream. In this sense, our theory is more general than the impulsively 
started scenario. 


5. Presence of swirl at pipe inlet 

Naturally, to maintain a uniform parabolic profile over the whole cross-section at 
the inlet is extremely difficult; the entry velocity is likely the combination of a 
normally parabolic distribution and a swirling component. We postulate that 


u = V = 0, w = w(r,0,t), and p = p{r,0,z,t) 


over time t G [0, e]. Here we have assumed that the 0-dependence of the pressure is 
a result of accelerating the flow into the inlet, possibly from a large fluid container, 
so that an equivalent body force F{r, 9, t) is generated. Let F^ and Fg be the 
components in the r and 9 directions respectively. Thus we must have dp/dr = F^, 
and dp/d9 = rFg. 

The equation of motion has a simple form of 


dw / d'^w 
dt dr'^ 


1 dw 



1 d'^w\ 

J.2 QQ2 } 


dp/p 
dz ’ 


w(r, 0, 0) = S'o(r,0), w{l,9,t) = 0. 


(5.1) 


The solution for the initial-boundary value problem can be written as 

^27r 

w{r^9^t) = / / ^0')K{r^0^r' ^0\t)dr'dO' 

Jo Jo 


/o ^0 

rt /*! n277 


0 ^0 ^0 


( dp/p 
V dz 


(5.2) 


{r' , 0', z' , T)K{r, r' , 0,0', t—t')dr'd9'dt', 
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where Green’s function K is given by 


1 ^ ^ 


^ (7o (AttiT) t/o ) 


m—1 


m—1 




(5.3) 


The constant crn,m is the mth positive zero of J„(cr) = 0 (Watson 1944). The 
derivatives of the Bessel function can be evaluated from the recurrence relations 

J^{a) = Jn-ii<j) - nJn{<j)/a and Jn+i(cr) = (2n)J„(cr)/CT - J„_i(cr). 

Suppose that the initial data near the pipe inlet can be expressed as 


S'o(r,6») = /3(1 -r 2 ) 7 ( 6 »). 


In the limit of vt^O and a constant axial pressure gradient, the leading order 
approximation by keeping the first sum in (lOl) is proportional to the angular 
integral of 7 ( 0 ). The resulting profile remains parabolic with modified strength. For 
the full effects of swirl, we must include the double sum and the analysis becomes 
much more involved in general though it is a matter of computation. 

The parabolic profile is not the only uni-directional time-insensitive flow achiev¬ 
able in a straight pipe. If the pressure gradient is allowed to vary as a non-linear 
function of x, a whole range of flow profiles, possibly with strong swirls, may be gen¬ 
erated (see, for example, Pfenninger 1961). The short-time evolution of these flows 
can be explained by our solution (03) and (03 as long as we have an adequate 
knowledge of the initial data wo(r) and the pressure near the pipe inlet. 


6. Numerical examples 


In this section, we present some numerical results for various initial data at the 
entrance of a long rigid circular pipe (the radius is assumed to be unity). The 
profile given in solution (03 is normalised by quantity —{dp/dz)/{4,p) = Pq: and 
we write the scaled velocity as 


w*{r,t) = w{r,t)/Po = wi{r,t) + WG{r,t), 


where the first term on the right stands for the effect of scaled inlet data, ic*, 
and the second term is induced by the Green function at constant axial pressure 
gradient. 


(a) Initial parabolic / boundary-layer profile 
In figures [T] and 01 we show two examples where the initial velocity distribution 


is a parabolic type or a Hagen-Poiseuille-Gouette type, while the pressure gradient 
is nominally constant to maintain a constant mass-flux through the pipe. If we 
assume that a motion is set up impulsively from rest, the flow development under 
the constant pressure gradient is in fact the solution of an impulsively started 
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motion from a parabolic velocity obtained by Szymanski (1932) where the scaled 
profile is given by (cf. equation (14.611 1 

WG = (1 - - 8^ exp ( - \lvt) ■ 

This is the solution presented in the middle rows in figures [T] and [21 The computed 
results show that initial data with parabolic distributions do not describe the flow 
regimes over the entry length. Specially, the viscous layers at pipe wall are mis¬ 
represented. In addition, parabolic profiles at pipe entry are difficult to generate 
in practice. If the entry profiles are of boundary layer type, as shown in figures [3] 
and |4l the developed velocity downstream of the entry resembles the well-known 
experimental measurements (see, for example, §2.5 of Tritton 1988). A parabolic 
velocity profile can only be realised in the limit oi vt ^ oo or at locations far 
downstream of the entry as long as the axial pressure variation is kept constant. 

(6) Entry profile with swirl 

We show one example of inlet flow with swirl with the natural period in the 
azimuthal direction: 


wl{r,e) = ^{l-ry/^(l + ^snfl{2ne)) 

As assumed in the present theory, the swirl is generated by an body force near the 
pipe entrance and the subsequent flow has a constant mass-flux or a constant axial 
pressure gradient. Because of the 9 rotational symmetry in Green’s function dO]), 
the second contribution in (lOl) is identical to the non-swirl cases just examined. 
Hence only the computational results of wj are presented in figure |5| where, clearly, 
the overall flow field has certain ’’symmetries” and is substantially more complicated 
to analyse. We shall carry out detailed investigation in separate studies. 

7. Conclusion 

We have shown that, by including du/dt term in the equations of motion, the uni¬ 
directional flow in a channel dSSl) or in a pipe (03) is space-time-dependent and 
exists over a short time when the equations are solved with non-zero initial data. 
The classical plane Poiseuille flow (12.31) or Hagen-Poiseuille profile (12.71) are best 
interpreted as an evolutional solution of an appropriate initial velocity in the limit 
of nt —>■ 0+. This point has been heuristically argued by Lam (2013) on the basis 
of exact vorticity solutions in subsequent development over longer periods 

must be governed by the full equations of motion where the non-linearity becomes 
dominant. It has long been known that, by neglecting the non-linear term (rt.V)rt, 
the approximate fluid motions are only valid for low Reynolds numbers, see, for 
example, §110 of Prandtl & Tietjens (1934). 

It is not difficult to understand why there do not exist absolute steady flows, 
independent of initial data, at arbitrary Reynolds numbers. Our flow regimes are 
clearly outside the physics of relativity, the concept of space-time entity enables us 
to recognise the conceptual inconsistency in the formulas (17:71) or (1731) . 
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Although the differences in the spatio-temporal solution (14.21) and the abso¬ 
lute steady flow (1^ appear to be insignificant with reference to experimental 
verifications, particularly for fluids of small viscosity, the hydrodynamic principles 
underlying the two fluid motions have important ramifications. For instance, the 
linear stability analysis for the non-dimensional parabolic profile shows that the 
basic flow is stable to all infinitesimal disturbances (see, for example, Lam 2014). 
By this counter-example, one casts serious doubts on the relevance of the linear 
stability theory in explaining the process of laminar-turbulent transition. Could it 
be something fundamental which may have gone wrong? If the basic flow evolves in 
space and time, what is the meaning of its stability or instability when the Reynolds 
number itself is varying? 


04 May 2015 
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Figure 1. Flow regimes near pipe entry by uni-directional flow approximation. Entry 
profile is supposed to be a parabolic entry velocity or a Poiseuille-Couette profile. 
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Figure 2. At large time, the velocity field driven by dpjdz becomes dominant and tends 
to a parabolic profile in the limit of vt —^ oo. In the limiting case, it is the Green function 
G(r, s, t) which determines the large-time decay behaviours. 
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Figure 3. Flows over the pipe entry length from two boundary-layer profiles. 
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Figure 4. The initial boundary layers decay in large time. 
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Figure 5. Velocity profiles at 8 angular locations are shown in plots (a) to (h) for vt = 0.01, 
0.02, 0.03, 0.05 (black lines) compared with parabolic profile and boundary layer without 
swirl (red lines). 
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